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^ ■ Abstract: Preliminary PHENIX data on Au + Au reactions at = 130 AGeV suggest j. 

that p yields may exceed tt" at high pT > 2 GeV/c. We propose that jet quenching in 

. . central collisions suppresses the hard PQCD component of the spectra in central + ^ o 

OO . reactions, thereby exposing a novel non-perturbative component of baryon dynamics. We C ^ 



suggest that baryon junctions provide a possible explanation of the anomalous component. 



^\ ' We predict that the p > n anomaly is limited to a finite window and decreases with 



increasing impact parameter. 



D I 1. Introduction 

Preliminary RHIC data on ^/s = 130 AGeV Au + Au reactions have revealed a number of 
^ I qualitatively new phenomena at moderate high pT ^ 2 — 6 GeV/c. 

' 1. The high pT spectra of vr'^ in the 10% central collisions were found by PHENIX to be 

suppressed by a factor 3-4 relative to PQCD predictions scaled by nuclear geometry 
(r^B(b), the number of binary collisions at impact parameter b). In contrast, the 
inclusive charged particle spectrum was reported to be suppressed only by a factor 
of 2 0. 

2. Non-central collisions STAR measured elliptic flow consistent with hydro predictions 
in the soft pi region. However, this asymmetry saturates at V2 ~ 0.15 above pT > 
2GeV/c |,|]. 

3. Even more surprisingly PHENIX reported 0] that the high > 2 GeV flavor com- 
position of positive and negative hadron yield may actually be dominated by anti- 
protons and protons respectively! 
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In addition to the results listed above, there are already global indications from pre- 
liminary STAR data ||5| of novel baryon production dynamics in ^ + ^. In central collisions 
the midrapdity p/p ~ 0.65 and dN^ /dy ~ dNP/dy - dN^/dy ~ 15. Thus a very signif- 
icant baryon number transport (stopping) five units of rapidity from the fragmentation 
regions and copious BB production seems to have been observed. This is consistent with 
expectations based on a non-perturbative baryon-junction mechanism for baryon number 
transport In this paper we propose that the baryon/meson anomaly together with the 
different degree of suppression of vr" and inclusive charged hadrons in central collisions may 
be an unexpected new consequence of jet quenching which unveils a novel non-perturbative 
component of baryon production dynamics 0]. 

Our approach is to extend the two component hybrid model introduced in Ref. Q 
considering an anomalously large contribution to p production at moderate p^ motivated 
by baryon junction dynamics 0. The baryon junction picture implies a simple relation 
between the mean inverse slopes of pions and protons, Tp ~ ^/2>T^^ (in the hassles limit), for 
the soft to moderate pT non-perturbative part of the spectrum. We compute jet quenching 
using the non-abelian energy loss from the Gyulassy-Levai- Vitev (GLV) formalism 
Our analysis suggests that the enhancement of p/tt" ratio is limited to a finite px range 
2-5 GeV and to central and semi-central collisions. Beyond this range, the hadron ratios 
are predicted to converge to the jet quenched PQCD base. 



2. Reference flavor composition in PQCD 

The standard PQCD approach expresses the differential hadron cross section in p -|- p — > 
h + X as a. convolution of the measured structure functions /(xq,, QDa/p for the interacting 
partons (a = a, 6), with the parameterized fragmentation function Df^i^{z^Ql) for the 
leading rescattered parton c into a hadron of flavor h and the elementary parton-parton 
cross sections da^"'^^^'^^ / di. The initial parton broadening can be accounted for by a 
normalized fcT-smearing distribution. A simple functional form is the Gaussian smearing 
/(/ct) = e~^T^/^^T:^ /'K{k'^) where (fc^) ~ 0.6 — 1.0 GeV^ in pp collisions. The invariant 
hadron inclusive cross section is then given by 

^h-^ J dZcdXadXb J d'^kTad'^kTbfO^Ta)fOiTb)fa/p{Xa,Ql)fb/p{Xb,Qb) 

abed 

^ ^^(ab^cd) 

Dh/cizc, Ql) — IT + n + t) , (2.1) 

vrz^ dt 

where Xa,Xb are the initial momentum fractions carried by the interacting partons, Zc = 
Ph/Pc is the momentum fraction carried by the observed hadron. 

Fig. la shows comparison between the PQCD calculation and the negative hadron 
multiplicities as measured by the UAl experiment at ^/s = 200 GeV. For our first fit 
we used the Gliick-Reya-Vogt structure functions (GRV 94) and Binnewies et al. (B) 
fragmentation functions. For our second fit we used the CETQ 5M structure functions and 



the Kniehl et al. (KKP) fragmentation functions |1C]. Eq.( |2.lD reproduces well the shape 
and norm of the data for px > 1 GeV. 
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Figure 1: a) Negative hadron cross sections are shown as a function of px- Experimental 
UAl data at ^/s = 200 GeV is compared to PQCD calculations with GRV94 pdf + B 
fragmentation functions (Q^ = v\l^) and CTEQ 5M pdf + KKP fragmentation function 
(Q2 ^ 2p|) where (fc|) = 0.8 GeV^ in both cases, b) The PQCD predictions for the 
p/tt+j p/tt^, 7r^/7r+ and p/p ratios as a function of are shown for -^/s — 130 GeV. 

Fig. lb, however, shows that although PQCD fragmentation can fit the 7r^/7r+ and even 
get quantitatively correct p/p ratios for px > 3 — 4 GeV the p/vr"^ and p/vr^ differ from 
the values suggested by data (~ 1) by an order of magnitude. Even large uncertainties 
in the gluon fragmentation functions (~ 3) Q cannot account for such a discrepancy. 
The observed large (anti)proton enhancement in Au + Au collisions in the region of px = 
2 — 4 GeV ^ signals of interesting interplay between competing possibly novel physical 
effects. We propose that the mechanism which accounts for the copious production of 
baryons and anti-baryons might be related to baryon junctions. Other possibilities include 
relativistic hydrodynamics and boosted thermal sources. 
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3. Medium effects and baryon junctions 

We first discuss the medium effects on the high px particle spectra. The medium induced 
gluon radiation spectrum is calculated as in and modifies the fragmentation function 
^hjJy^c) Eg. (p. ID as follows 

'1 



dx P[x) 



1 



1 



-D 



X 



h/c 



1 



X 



X 



h/g 



(3.1) 



where x = cv/Ejet, P{x) is the gluon radiative probability, and the second term in Eq.(^ 
accounts for the gluon feedback in the system. Transverse expansion was shown not to 
affect the azimuthally averaged quenching pattern |^] 

The string and baryon junction description of the soft non-perturbative part of the 
spectrum is parametized in terms of the mean inverse slopes of hadrons Tq as in The 
baryon junction picture, however, suggest different Tq and Tq. Large Tq may arise from the 
small ~ 1/2Tq ~ 1/4 fm) intrinsic (non-perturbative) spatial structure of Junction- 
anti-Junction loops and the smaller junction trajectory slope a'j « ot'^/'i The latter 
implies that the effective string tension is three times higher than l/(27ra^) ~ 1 GeV/fm 
leading in the massless limit to 



'T/J 



3(p2 



rpp 

^0 



(3.2) 
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Figure 2: a) Inclusive charged hadrons and b) neutral pions are plotted versus pT nor- 
malized by the corresponding scaled PQCD multiplicities. Curves reflecting the effect of 
non-abelian energy loss driven by dN^/dy = 500, 1000, 2500 are shown. 

Eq.(|3.2|) helps reduce the parameters that enter the phenomenological soft component to 2: 
Tq ~ 210 — 220 MeV and dN~^/dy ~ 600 in central collisions at midrapidity. Assuming that 
similar string and baryon junction dynamics drives pp and AA collisions, we solve for the 
remaining parameters from the the existing data |jll| (p/vr^ ~ 7 — 8%, /tt^ ~ 12 — 15%) 
or from charge and strangeness conservation (vr-l- /vr^ ~ 1, /K~ ~ 1 — 1.1) and baryon 
transport @ {dN^/dy ~ /3Zcosh(l - aij(0)y)/sinh(l - ai3(0)ymax) 10 - 15, given 
l^nax = 5.4, /3 ~ 1). 
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4. Result and conclusions 

In Figs. 2a and 2b we compute the quenching for inclusive charged hadrons and neutral 
pions. The ~ 2 difference in the suppression factors in the = 1 — 4 GeV window is largely 
due to the non-perturbative baryon component of dN^h- The comparison of the quenching 
factors to data ^ suggest initial gluon rapidity density dN^/dy ~ 500 — 1000. Fig. 3 
presents the central result of this paper. The anomalous baryon/meson ratio exceeds unity 
in the pT = 2 — 5 GeV window in central collisions. It is here interpreted as a combined 
effect of non-perturbative baryon junction dynamics and jet quenching that suppresses the 
PQCD component of pions. Fig. 3 suggests that this result may be more readily observable 
in p/tt~^ than in p/'K~ . The uncertainties in the mean inverse slope of pions and protons 
as well as the degree of quenching may affect Rb max by as much as 30 — 50% (also seen in 
Fig. 3), without however qualitatively changing our predictions. Boosted thermal sources 
calculation with vt = tanhr/^ ~ 0.6 and Tj ~ 160 MeV 



dN, 



dyd'^p^ 



TTT-T Ki 



TTT-T cosh T]^ 



Pt sinh r/^ 



(4.1) 



may produce similar low p^ behavior but Rb is a monotone function and saturates at 
Rb = 2 from spin counting. We predict that at high > 5 — 6 GeV Rb decreases again 
below unity to its PQCD dominated base. We also predict that in going from central to 
peripheral collisions the baryon/meson anomaly disappears since there is no jet quench- 
ing to suppress the perturbative pions and the Rb ratio resembles the pp measurement. 
Preliminary PHENIX data seems to support the suggested centrality dependence. 
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Figure 3: The p/vr^ and p/7r+ ratios are plotted versus px for 4 different impact parameters 
b = 0, 4, 8, 12 fm. The top and bottom row of figures illustrate plausible uncertainties in 
the mean inverse slope of tt and p. The ratio of our fits to preliminary PHENIX data as 
well as boosted thermal source calculation are shown for comparison. 
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